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We use the VESPA algorithm and spectra from the Sloan Digital Sky Survey to investigate the star formation history of the host 
galaxies of 257 Type la supemovae. We find Scr evidence for a short-lived population of progenitors with lifetimes of less than 180 
Myr, indicating a Type la supernova channel arising from stars in the mass range ~3.5-8 Mq. As standardizeable candles. Type la 
supemovae play an important role in determining the expansion history of the Universe, but to be useful for future cosmological 
surveys, the peak luminosity needs to be free of uncorrected systematic effects at the level of 1-2%. If the different progenitor routes 
lead to supemovae with even moderately small differences in properties, then these need to be corrected for separately, or they could 
lead to a systematic bias in future supemovae surveys, as the prompt route is likely to increase in importance at high redshift. VESPA 
analysis of hosts could be a valuable tool in this, by identifying which progenitor route is most likely. 

Key words, supemovae: general — galaxies: stellar content — distance scale 



1. Introduction 

The relationship between apparent peak brightness and redshift 
of Type la supemovae (SNe la) depends on the cosmological 
model; this has provided the most direct evidence for the accel- 
erated expansion of the Universe ( |Riess et aL] |1998[|Perlmutter| 
[et all [1999 ). Current SN la surveys such as SNLS (|Astier et al. 



2006 ), ESSENCE ( |Wood-Vasey et al.||2007) , and GOODS-SN 
( |Riess et al. ,2007[ ) are contributing to current constraints on cos- 
mological parameters, and SNe la will continue to be important 
for cosmological constraints in the next generation of surveys 
such as the J DEM candi dates ADEPT, DESTI NY ( .Bendford &, 
[Lauer||2006| ) and SNAP ( |Aldering etaT||2004| . 

Type la supemovae are interpreted as the thermonuclear ex- 
plosion of a white dwarf that has reached the Chandrasekhar 
mass and thus has become unstable, probably through accretion 
from a companion star (the single-degenerate scenario) or merg- 
ing with another white dwarf (the double-degenerate scenario). 
However, no fully consistent model of a SN la explosion has yet 
been built. 

The natural scatter in SNe la peak luminosities covers 
roughly one magnitude; the rms peak luminosity is 0.45 mag 
after excluding outliers. Empirical correlations based on light 



curve shape ^ Phillips 1993) or intrinsic color ( [van den Bergh 



1995[ [Tripp & Branchy |1999j) allow r eduction of the intrinsic 



scatter to about 0.13 mag (Astier et al. 2006), making them us 



able for cosmological measurements. 

However, it is not yet known how much of the residual scat- 
ter is correlated with physical parameters that could evolve with 
redshift, and thus bias the measurement of cosmological param- 
eters. In order for SNe to be useful for constraining dark en- 
ergy at the level expected in future SN satellite experiments, 
the evolution of luminosity at a given light-curve shape over 



the probed redshift range must be less than 1-2% ( [Howell et al.| 
[2007; Sarkaret al.j ,2008) . 

There are several observational indications for a range of de- 
lay times between the birth of the progenitor system and the ex- 
plosion of the SN la, leading some authors to envision the ex- 
istence of at least two different populations of SN la populating 
slightly different regions of stretch-brightness parameter space 
(Sullivan et al. 2006). It is not yet known if they are described 
(at the percent level) by the same Phillips relation; if they corre- 
spond to different production channels there is no reason for the 
Phillips relations to be identical. In any case, the SN population 
may not be a one-dimensional family, but may depend on other 
parameters (at an as-yet unknown level) such as metallicity and 
delay time. This would result in a source of scatter in the SN 
la Hubble diagram that could be reduced with measurements of 
those additional parameters. Moreover, if the Phillips relations 
are different and the relative numbers of SNe in the two putative 
populations evolve with redshift, the values of cosmological pa- 
rameters derived from the Hubble relation will be biased if we 
do not determine population-dependent corrections. More gener- 
ally, any extra parameter (age, metallicity, production channel...) 
in the Phillips relations that is correlated with redshift will bias 
cosmological measurements. It is therefore essential to assess 
the size of such effects. Measuring them by looking for congela- 
tions between Hubble diagram residuals and intrinsic properties 
such as delay time, progenitor metallicity, etc. could help correct 
for them (at least statistically if not event-by-event) in future sur- 
veys. 

The brightest supernova events only occur in actively star- 
forming galaxies (Hamuy et al. 1995, 1996a,b), suggesting 
prompt explosions, while under-luminous events are most often 
found in spirals and E/SO galaxies, whose old stellar populations 
would suggest delayed explosions ( [Howell[[2001| ). Mannucci et 
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al. (2006) and |Scannapieco & Bildsten| ( |2005| l have pr oposed a 
two-component model for SNe la, and several authors (Sullivan' 
et al. 2006 Mannucci et al. 2005 |l have shown that the super- 
nova rate can be expressed as a sum of a term proportional to the 
total mass of the galaxy and a term proportional to the recent star 
formation rate. In the Mannucci et al. model, some of the super- 
novae would explode several Gyr after the birth of the progenitor 
system, while others would after a fraction of a Gyr Moreover, 
those two populations have different luminosities, the "prompt" 
component being brighter with broader light-curves (Sullivan et 
al. 2006). The prompt component will dominate at higher red- 
shifts when the Universe's age (or the time since star formation 
began) was less than the lifetime of the longer-lived progenitors. 
Determining the relative numbers of supernovae in the two pop- 
ulations is the first step in understanding any possible bias these 
populations might cause. Given the 13% scatter in the calibrated 
peak luminosities of supernovae, measuring systematic effects 
in the Phillips relation at the percent level will take samples of a 



few hundred supernovae ( Sarkar et al. 2008 1 



Hamuy et al. (2000) used 62 SN la host galaxies to study the 
impact of host morphology, magnitude and colors on the decline 
rate Amis, which allows one to estimate the SN peak luminosity. 
They found a correlation with age but not with metallicity (see 
also the erratum to their paper). However, their sample was very 
small and most of their estimates of age and metallicity were 
based on photometry only, without spectra of the host galax- 
ies, and therefore their accuracy was limited. Moreover, although 
they investigated various environmental effects, their methodol- 
ogy was not sensitive to a second parameter in the Phillips re- 
lation, since they used the decline rate as a "reddening-free and 
distance-free estimate of the SN peak brightness", and thus as- 
sumed a priori the universality of the relation. 

Gallagher et al. (2005) carried out a similar analysis with 
spectra of 57 SN la hosts, and put tentative constraints on the 
SN progenitor lifetime using an estimate of current-to-average 
star formation rate. They claimed to see hints of both a bimodal 
behavior and a lower limit of the progenitor lifetime. They ad- 
mitted that their findings were rather inconclusive. 

[Sullivan et aL| (2006) used 100 SNe from the SNLS and 
broadband spectral energy distributions of the host galaxies to 
estimate stellar masses and star formation rates. They found a 
component proportional to the stellar mass, and a component 
proportional to the recent star formation rate, averaged over the 
last 0.5 Gyr. 

Our study improves on these earlier papers by using a larger 
sample of SNe, with spectra of their host galaxies from the 
Sloan Digital Sky Survey (SDSS; York et al. 2000). We use a 
sophisticated stellar population code called VESPA ( [Tojeiro et 



[aLl[2007| ) which allows us to determine the stellar formation his- 
tory of the hosts. We also determine the star formation history 
of a large sample of normal galaxies from the SDSS as a con- 
trol. Differences in those star formation histories will yield in- 
formation on SN progenitor lifetime: short lifetimes for instance 
should statistically enrich the host sample in galaxies with recent 
star formation. 



2. Host galaxies and reference sample 

We gathered a sample of about 1300 confirmed SNe la from lAU 
circular^ the CfA supernovae lisj^and the SDSS-SN public list 

' http : //www ■ cf a . harvard . edu/iau/cbat . html 
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Fig. 1. Distribution of the number of stellar populations recov- 
ered by VESPA from the host sample (gray) and the control sam- 
ple (white). The two distributions are very similar (see text). 



of supernova^ We cross-re ferenced this fist with the SDSS DR5 
( Adelman-McCarthy et al. 2007 1 spectroscopic survey of galax- 
ies: 256 galaxies with spectra were identified as SN la hosts, 
corresponding to 257 supernovae (one galaxy hosted two super- 
novae). The list of hosts used in this paper is available onlin^ 
A large fraction of our sample comes from surveys like SDSS 
(104 SNe), LOSS and LOTOSS (49 SNe for the two surveys). 
The rest have various origins (Pucketf]: 17 SNe, Pollas (1994) : 
12 SNe, etc.). 

The detection efficiency of this sample is unknown, as it de- 
pends in detail on the way in which the SNe were found. To 
account for the selection function of SN la discovery, we also 
process a control sample of 10^ DR5 galaxy spectra, chosen 
randomly in the survey release, and weighted to reproduce the 
redshift distribution of the host sample — this is the param- 
eter which could most significantly bias delay time measure- 
ments. We plan to handle these effects more accurately with a 
full Monte Carlo treatment in a future paper. 

As a cross check, we did the same analysis, keeping only 
supernovae at z < 0.1. The reduction in the sample size (from 
257 to 190) increases the error, but does not change the result 
(see table 2 below). Other effects will be discussed in section 4. 



3. Reconstructing the star formation and metallicity 
history of SN host galaxies 

The spectrum of a galaxy is a superposition of spectra of sin- 
gle stellar populations which formed at a given age with a given 
metallicity. Since it is not possible to recover the star formation 
and metallicity history with infinite precision (e.g. Jimenez et al. 
|2004[ ), it is only sensible to attempt to recover the star forma- 
tion and metallicity history with a certain time resolution. The 
VESPA algorithm (Tojeiro et al. 2007| l does this, providing a 



|http : // cf a-www . harvard . edu/ iau/li sts/Supernovae . html ^http : //www . cometwatch . com| 
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detailed history only where the data warrant it. Note that broad- 
band colors are not sufficient to determine the star formation his- 
tories of galaxies, as they suffer from significant age-metallicity 
degeneracies (e.g. Jim enez et al.||20(j4] l. 

In brief, VESPA uses singular value decomposition to cal- 
culate the number of significant components in the spectrum of 
a given galaxy. VESPA then uses an algorithm to determine the 
best-fitting non-negative values of the star formation fractions. 
Extensive tests of the performance of VESPA on synthetic spec- 
tra as a function of wavelength coverage and signal-to-noise ra- 
tio can be found in Tojeiro et al. (2007 1, along with a study of 



age-metallicity degeneracy. To limit the search to a manageable 
amount of parameters, and because currently available spectra 
never have the quality or spectral range to justify going be- 
yond this choice, VESPA's finest resolution consists of 16 age 
bins, logarithmically spaced in lookback time between 0.002 
and 14 Gyrs. Specifically, the lower limit in age of the 16 bins 
are: 0.002, 0.02, 0.03, 0.0462, 0.074, 0.115, 0.177, 0.275, 0.425, 
0.6347, 1.02, 1.57, 2.44, 3.78, 5.84, and 9.04 Gyr 

VESPA chooses the number of stellar populations to model 
depending on the quality of the data. The SDSS galaxy spectra 
typically allow between 7 and 10 age bins in both the SN host 
sample and the control sample, showing there are no significant 
differences in the two samples (Fig. [Tj see also Tojeiro et al. 
2007), although there is non-zero star formation in only 3-5 of 
those bins. The relatively large number of bins implied by Fig.[T] 
is an indication that in most cases the star formation which is 
recovered is localised to narrow time intervals. The metallicity 
for each population is a free parameter, so there are as many 
metallicity values recovered as there are star formation fractions. 

Following Sullivan et al. (2006), we parametrize the total 
supernova rate with a two-component model, SNR - aMt + 
fiMrecent, whcrc M» is the total stellar mass and Mrecent is the 
mass formed in some time range that we can vary, a and fi re- 
flect, respectively, the SNR per unit stellar mass of an old pop- 
ulation of progenitors (roughly proportional to the total stellar 
mass), and the SNR per unit stellar mass of a young popula- 
tion of progenitors (proportional to the mass in recently formed 
stars). Because we do not have a calibration of the selection func- 
tion of the supernova surveys we have used, we normalize with 
our control sample, representative of the whole SDSS spectro- 
scopic catalog (our reference sample), and can determine SNR 
only up to a proportionality constant. That is, we fit the data to a 
model 



SNR 

~m7 



c 1 



a Mt 



(1) 



where C is an unnormalized proportionality constant. Thus the 
quantity /3/a quantifies the fraction of the supernova rate with 
progenitor stars that formed in the chosen time range. The con- 
stant C is the product of the "slow" rate A introduced in the 
Mannucci et al. two-component model, and an unknown global 
SN detection efficiency factor averaged over all SN surveys 
which contribute to our reference sample. 

We estimate the SNR per unit stellar mass, SNR/Mt as a 
function of as follows. For any unbiased sample of galax- 
ies, the observed SNR per unit stellar mass is equal to the num- 
ber of supernovae divided by the total stellar mass in the sample, 
multiplied by an efficiency factor e: - f^g^- We can di- 
vide them into bins of their star formation fraction in the recent 
time range, , as determined by the VESPA analysis. Within 

each bin we can therefore calculate by simply counting the 
number of SNe in each subset and dividing it by the total stellar 



Age (Myr) 




Error on yS/ff 


Significance (cr) 


0-180 


454 


78 


5.8 


180-660 


56 


16 


34 


660-2440 


18.4 


3.5 


5.2 


2440-13700 


-3 


1.0 





Table 1. Fits of supernova rate per unit stellar mass oc 1 -n-^ 
for different ranges of star formation lookback time for the /3 
component. 



Age (Myr) 


/3/a 


Error on/} /a 


Significance (cr) 


0-180 


347 


75 


4.6 


180-660 


43 


20 


2.2 


660-2440 


46 


8 


5.8 


2440-13700 


-3 


1.0 





Table 2. Same as table 1, considering only SNe with z < 0.1. 



mass in galaxies with this value of — The latter is calculated 
from the much larger control sample. 

We repeated this exercise for a variety of ages for the young 
component, and also varying the dust model details in VESPA 
and the bin boundaries. We find a significant contribution to the 
SN rate from recent stellar populations (in the sense that /3/a is 
significantly positive), where we vary the definition of "recent" 
between 74 and 180 Myr. The contribution of stars in the 180- 
250 Myr range is lower by at least a factor of five. The value of 
/3/a is robust to setting the boundary anywhere between 74 and 
180 Myr, and thus here we quote the most conservative value 
of 180 Myr, since starburst duration and internal degeneracies in 
VESPA could make age bins "leak" one into another 

We illustrate this in Figure [2] We have divided our 16 time 
bins into four broader bins, and asked for the correlation of the 
supernova rate with the fraction of star formation that occurred 
in each bin. The four panels show the un-normalized Type la 
supernova rate per stellar mass as a function of the fraction of 
stellar mass formed each of these broad time bins [0-180] Myr 
(top-left panel), [180-660]Myr (bottom-left panel), [0.66-2.44] 
Gyr (top-right panel), and [2.44-13] Gyr (bottom-right panel). 
The resulting best-fit value of /3/a and the corresponding corre- 
lations are given in Table 1 . 

For the youngest range, 0-180 Myr, we measure /3/a - 454 ± 
78, which is five-sigma evidence for a short duration component. 
Since /3/a in older bins is significantly lower, the positive signal 
we find in the most recent age range cannot be due to correlation 
("leakage") from older bins. 

Our results should be robust against possible spectroscopic 
calibration errors: because we compare the host population 
to a control sample of spectra taken with the same telescope 
and instrument, and processed with the same spectroscopic 
pipeline and the same star formation history recovery algorithm 
(VESPA), systematic errors anywhere in the chain would be 
shared by the two samples. 

The selection function of our SN sample is also irrelevant as 
long as it is not strongly dependent on host properties which are 
correlated with SFR; the effect of possible efficiency biases in 
the SN host sample will be discussed in more detail in §4 below. 
We do not consider supernovae without a SDSS host but this is 
fully consistent with our approach of comparing to a SDSS ref- 
erence population : our reference sample is, by definition, com- 
plete, and we compare, among that sample, hosts without SN to 
hosts with SN. 
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Fig. 2. Type la supernova rate per stellar mass, unnormalized, vs. fraction of stellar mass formed in the four time intervals indicated 
in each panel. The dashed line is a fit to a dual component model S NR = axMt+/3x Mtime range, for each time range. In the absence 
of a stellar population of a given age contributing to the supernova rate, the curve would be consistent with flat (i.e,, /3 consistent 
with zero). We find that /3/a is non zero at the five-sigma level for the time intervals [0-180] Myr. Shaded areas are ±2cr fit values. 



A bias will be introduced in our results if there is a cor- 
relation between galaxy mass and specific star formation rate. 
Indeed, most massive galaxies in the present-day universe have 
little star formation; most of the star formation in the present-day 
universe is in low-mass galaxies (Heavens et al. 2004). We have 
tested this bias with a Monte-Carlo simulation, as follows. 

We generated an artificial SN sample in our control sample 
following the rule a x M» + f5 x Miecent- In the absence of cor- 
relations between recent star formation and mass in this refer- 
ence sample, applying our method recovers the input /3/a values. 
However, our reference sample has an anti-correlation between 
mass and recent star formation. With the input /3 set equal to 
zero, we recover jS/ff = -3+2, compatible with zero (dashed line 
in Fig. [3]l. Using /i/a - 700 (dotted line in Fig. [3]l, we recover 
/i/a - 465, 2/3 of the simulated value (and close to our observed 
value). Therefore we conclude that the anti-correlation between 
galaxy mass and star formation rate causes us to underestimate 
the contribution of a prompt component to the supernova rate. 

The ratio /3/a is compatible with previous estimates, for 
which the "recent" SFR is estimated in general from colors, 
broadband SED fitting, core-collapse SN rate or cosmic SFR. 
These results represent an average over half a gigayear, and are 
thus only a rough match to our results. We can roughly con- 
vert our mass estimate Mi so {M recent for the last 180 Myr) to 
recent SFR through Migo = 180 x 10* j^|^Mq. With this, the 

Neill et al. (2006) values (a "slow" rate A = 1.4 + 1.0 x lO"''* 



SN M-i yr-i and a "prompt" rate B = 8.0 + 2.6 x lO""* SN 
(MQyr-i)-iyr"') yield /3/a ^ 300, the Sullivan et al. (2006) 
values (5.3 + 1.0 x lO^'^ SN M^' yr"' and 3.9 x lO""* SN 
(M0yr"')"'yr"') yield /3/a ^ 40, and the two values quoted by 
Scannapieco & Bildsten (2005) yield /3/a ^ 300 and 150, re- 
spectively. However since "recent" star formation is estimated 
differently in each case, and on an difference recent time range, 
the figures are not expected to match closely. 

4. Discussion and conclusion 

Our result would be sensitive to any systematic effect enrich- 
ing the SN host sample in blue galaxies (i.e., those with large 
Mrecent/Mt), foi rcasons unrelated to SN physics (bias in ef- 
ficiency, or bias in the monitored galaxy sample for targeted 
searches). Some targeted supernova searches have focused on 
star-forming galaxies (e.g., Richmond, Filippenko, & Galisky 
1998). In our sample, the main targeted search is the Lick 
Observatory Supernova Search (LOSS, Filippenko et al. 2001). 
There is no indication this search is biased in this way, and our 
results change insignificantly if we remove those hosts. 

Blue galaxies are fainter, and one could naively expect to de- 
tect SNe more easily in those faint hosts, but with modern image 
differencing techniques, the effect should be minor. One also ex- 
pects such a bias to occur in surveys that rely on spectroscopic 
typing of SNe, since the fainter, low-stretch SN, that occur pref- 
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Recent star formation fraction 

Fig. 3. Estimation of the possible bias that can be introduced by 
star formation today taking place preferentially in smaller mass 
galaxies. Using our SDSS control sample we have generated su- 
pernova using jS/ff = and /3/a = 700. The recovered lines 
show that for the case with /3/a - (dotted line) we indeed find 
no signature of supernova, while for the case /3/a - 700 we do 
recover supernova activity at the 2/3 level (dotted line) of the 
input value (J3/a ^ 465), thus our conclusions are conservative. 
See more details in text. 



erencially in bright galaxies, will be more difficult to charac- 
terize spectroscopically. Star-forming galaxies tend to produce 
brighter events. On the other hand, star-forming hosts tend to 
be dustier, making supernovae harder to detect or characterize. 
If such an efficiency effect had an impact on our analysis, in the 
sense of enhancing or simulating what we observe, we would ex- 
pect the host galaxy sample to exhibit lower typical masses than 
the reference sample, or lower dust content. VESPA yields an es- 
timate of dust content and luminosity of the hosts. We see no sig- 
nificant difference in luminosity and dust distributions between 
the host and control samples. Rather, the host stellar mass dis- 
tribution is shifted towards slightly higher masses, as we would 
expect given that the SN rate should increase with the mass. We 
are thus confident this effect does not significantly bias our re- 
sults. 

Such an efficiency effect would also be more important at 
high redshift. As we noted above, restricting our sample to SNe 
occurring at z < 0. 1 does not change our findings. In our next pa- 
per, we plan to use the lightcurve characteristics in our analyses 
for those objects for which they are available. 

Following Mannucci et al. (2005, 2006) and Sullivan et al. 
(2006), we have shown that SNe la can occur through short- 
lived progenitors, hinting at a variety of stellar evolution paths 
with different lifetimes. We have given the first estimate of the 
lifetime of the "fast" component, by reconstructing the star for- 
mation history of SN host galaxies and finding an increased con- 
tribution to the SN la rate from stars evolving in less than 180 
Myr. 

Such a short time delay strongly constrains the nature of pos- 
sible progenitors. They must be stars that evolve fast enough, i.e. 
with a mass above ~ 3.5Mo, but must be below the super- AGB 
mass threshold (about 8 Mq) above which one gets electron- 

2008] l. 



capture supernovae (Poelarends et al 



Pinsonneault & 



|Stanek| ( 20061 have also suggested that a significant fraction 



of binaries are twins (i.e., pairs of stars with essentially iden- 
tical masses), and that such twin binaries could produce a short 
(< 0.1 Gyr) path to SN la. Considering common envelope evolu- 



tion phenomena, [Pinsonneault & Stanek| ( |2006| l argue that such 
twin systems could yield double degenerate SNe la in a way that 
would be both fast and efficient (see also H achisu et al.|2007] l. 

Are there enough high-mass progenitors to account for the 
observed SN la rate? These progenitors have to have masses be- 
tween 3.5M0 (in order to explode within 180 Myr) and EMq. 
Only a fraction of these stars ^ will actually explode as a SN 
la progenitor. We take into account five factors: the fraction 
of stars in binaries (/„), the fraction of the binaries both of 
whose components lie in the range 3.5 to 8 M© { ft,), the frac- 
tion of stars at a suitable separation for mass transfer (/c), the 
fact that every binary yields a single explosion, (/^), and an 
overall efficiency (77^, as not all possible progenitors may ex- 
plode). Maoz (2007) has estimated the first four factors, and 
finds /„ e [2/3,l],/i e [1/6,1/3],/, e [1/4,1/2] and./, = 1/2. 
Crudely multiplying these factors together gives the fraction of 
objects in the appropriate mass range that explode as prompt 
SNe la:/^ e [0.014, 0.083]/;/,. 

The fraction of stars which will explode as a SN la pro- 



genitor, /g is then given by 



A'3.5-8 



where A^^ is the total 



number of SNe la from the fast route, and N^.s-^ is the total 
number of stars in the correct mass range. Our result only al- 
lows us to estimate the ratio of the two components. We can 
however use published values of A in order to infer an abso- 
lute value of B from our ratio — we call this value B'. Using 
A = 1.4 + 1.0 X IQ-^'^SNyr^M^^ as published by Neill et al. 
(2006) gives B' = 6.3 ± 4.7 x lO^i^^A^yr-i M^'. We can now 
estimate A^B - B' x 180MyrsxMi8o andA^3.5_g - 0.0157*Mi8o, 
assuming a Salpeter IMF of the form N{m) cx m^^-^^. This gives 
/g = 0.073 + 0.053, within the values given by Maoz (2007). 
Using a value of A = 5.3 ± 1.1 x lO-'"* from Sullivan et al. (2006) 
gives a value of B' = 2.41 ± 0.65 x 10"" and //j = 0.28 + 0.07. 
Similarly, A - 4.4;^[-^ x 10"'"* as published by Scannapieco & 
Bildsten (2005) results in roughly B' = 1.99 ± 0.69 x 10 " and 
fl3 = 0.23 ± 0.08. 

With the slow A rate of Neill et al (2006) there is complete 
consistency with the theoretical expectations of SN la rates from 
Maoz (2007). However, with the higher A rates of Sullivan et al 
(2006) or Scannapieco & Bildsten (2005), there is some tension 
with our results, which would require high efficiency. It is hard 
to know how worried one should be about this: firstly we are 
relying on external measurements of the delayed rate, and sec- 
ondly an excess of SNIa explosions with respect to the predicted 
number of progenitors is observed in a large number of SNIa 
studies (Maoz 2007). This tension can be alleviated in a num- 
ber of ways — some of which discussed in the above paper — but 
generally indicates that the efficiency of the mechanism which 
produces SNIa explosions must be high. 

A crucial question for the use of SNe la as standard can- 
dles in cosmology is whether these different routes yield objects 
which are standardizeable to high accuracy via the same empir- 
ical corrections. Current data find no evidence for a difference 



Pamuy et al.][2000t[Sunivan et al.|[2003j|Bronder et al.|[2008) , 
but the requirements for using SNe la as a Dark Energy probe 
are stringent, and it will be important to establish this point ac- 
curately. VESPA should be able to assist directly with the correc- 
tion, first by determining the standardization for each route at the 
required percent level accuracy, and then allowing to apply the 
correct standardization at least statistically if not for individual 
supernovae. Furthermore, VESPA provides metallicity estimates 
for the hosts, which may correlate with peak luminosity, and thus 
allow us to reduce the scatter in the distance indicator. This will 
be the subject of a future study. 
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We plan to expand our sample by obtaining spectra of a 
larger number of SN hosts, allowing us to deconvolve the delay 
time function. Future papers will address more quantitatively the 
long duration component, the downsizing bias, the metallicity 
effects (see also Prieto et al. 2007), and stellar evolution models 
compatible with our findings. 
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